QUADRATIC HARNESSES, g-COMMUTATIONS, AND 
ORTHOGONAL MARTINGALE POLYNOMIALS 



WLODZIMIERZ BRYC, WOJCIECH MATYSIAK, AND JACEK WESOLOWSKI 

Abstract. We introduce the quadratic harness condition and show that in- 
tegrable quadratic harnesses have orthogonal martingale polynomials with a 
three step recurrence that satisfies a g-commutation relation. This implies 
that quadratic harnesses are essentially determined uniquely by five numerical 
constants. Explicit recurrences for the orthogonal martingale polynomials are 
derived in several cases of interest. 



1. Introduction 

Hammersley introduced harnesses on R" as the probabihstic models of long- 
range misorientation in the crystalline structure of metals. Several authors studied 
mathematical aspects of the concept: Mansuy and Yor analyze harnesses on 
R+, Williams analyzes harnesses with arbitrary discrete index set; see also |14| . 

isni, m- 

The class of random fields on M+ which we call the quadratic harnesses is re- 
lated to Hammersley's harnesses by mimicking the relation between the martingale 
and the quadratic martingale conditions. Such processes have already been studied 
implicitly by several authors, see the paragraph following Definition l2.2l In particu- 
lar, ref. gives a construction of the three parameter family of Markov processes 
with the quadratic harness property. Some of these processes are related to the free 
Levy processes, and some correspond to the non-commutative g-Gaussian processes 
introduced by Frisch and Bourret |17j and studied in [Sj. 

In full generality, the constructions of the quadratic harnesses on R_|_ are not 
yet well understood. In this paper we concentrate on uniqueness without dealing 
with the issue of existence. We show that quadratic harnesses are described by five 
numerical constants, which, under appropriate integrability conditions, determine 
the process uniquely. We study the related integrability properties of a slightly 
wider class of processes, improving earlier results of that type [HI Corollary 4], |10l 
Theorem 2], |32[ Theorem 2(1°)]. We show that martingale polynomials associated 
with a quadratic harness with finite moments of all orders lead to a g-commutation 
equation 

(1.1) [x,y], = I + Tx2 + (Ty2 + 6ix + ?7y. 
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where X, y are non-commutative variables, 

[x,y]9 -xy-qyx, 

and q,r],9,a,T are the numerical constants that describe the quadratic harness. 
This hints on more connections with non-commutative probability, g-Fock space 
constructions, and classical versions of non-commutative processes, see In this 
paper we use to prove that quadratic harnesses with finite moments of all 

orders have orthogonal martingale polynomials, and to derive their three step re- 
currences. The explicit three terms recurrences that are associated with quadratic 
harnesses include a four-parameter family of the polynomials in Section IT^ and a 
four-parameter family of q-orthogonal polynomials in Section 14.31 Based on one of 
the explicit recurrences from this paper, in |^ we extend the construction from 
to the construction of the quadratic harness which we call the bi-Poisson processes. 

Some of our arguments rely on the symmetries of the problem, but perhaps we 
did not explore the symmetries deeply enough. The family of quadratic harnesses 
that we study is invariant under the action of the translations and reflections of K, 
and the associated afline Hecke algebra is known to be associated with the Askey- 
Wilson polynomials, see (22]. Since the explicit recurrences that we found in Section 
^give sub-classes of the reparametrized Askey- Wilson polynomials, it is plausible 
that the theory of affine root systems 5^ might lead to additional progress. 

The paper is organized as follows. In Section [21 we define quadratic harnesses 
and state the main results. Proofs of the main results are given in Section |31 In 
Section 21 we derive explicit recurrences for the orthogonal martingale polynomials. 
We also introduce an operator technique that resembles umbral calculus as the tool 
that simplifies the proof of the quadratic harness property. 

2. Definitions and main results 

2.1. Harnesses and quadratic harnesses. Let {Xt)t>o be a separable square 
integrable stochastic process with cr-fields: Tg.u — (^{Xt : t E (0, s] U [u, oo)}. 
Consider the following two functionals 

(2.1) At,u^AtAX) = 



(2.2) At,„ = At,„(X) = 



u — t 
uXt - tXu 



u — t 

Both expressions appear in the statement of Theorem 12.21 below and the second 
one is the time-inverse of the first one. Namely, if {Xt)t>o denotes the time-inverse 
process {tXi/t)tyQ, then 

At,u{X) = Ai/„^i/t(X), Q<t<u. 

Definition 2.1 (|23). An integrable process {Xt)t>Q is a (simple) harness if 

E[A,,t|.F,,„] = A,,„ 

for every r < s < t < u. 

A trivial example of a harness is a Gaussian process with covariance (|2.5(l ; addi- 
tional examples follow Definition 12. 21 

As pointed out in |211, the harness condition is equivalent to the linearity of 
regression property 

(2.3) E[Xt|.F«,„] = at,s,uXs + bt,s,uXu, 
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where the coeificients are given by 

, , u — t t — s 

(2.4) at,s,u = , Ot.s,« = ■ 

u — s ' u — s 

It is clear that at,s,u+bt s u = 1 f^nd sat g y^-\-uht^s,u = t] these two identities will often 
be used. Since condition (|2.3|l is invariant under the time-inversion (Xt) i— > (tXi/f), 
another equivalent condition for a harness is 



E 



— Ar 0,, r < s < t < u. 



The general form of the covariance of a square-integrable harness is as follows. 

Proposition 2.1. If (Xt) is a square-integrable centered harness on M_|_, then there 
are constants cq, Ci, C2, C3 such that for s < t we have 

(2.5) E{XtXs)=ca + cis + C2t + C3st. 

Proof Multiplying H2.8|l by Xg and averaging we get E{XsXt) = at.s,uE,{Xf) + 
bt,B,uE{XsXu). Thus E{XsXt) = a{s) + (3{s)t is linear in t and E{XtXu) = 
at,s,uE{XsXu) + bt,s.uE{X^) = -f{u) + S{u)t is linear in t. 
This gives a functional equation 

a{s) + P{s)t = -f{t) + S{t)s 

which, by taking s = 0, gives j{t) ^a{0) + /3{0)t. Thus E{XsXt) = cq + ci.s + C2t + 
c^st. □ 

For at centered standardized square-integrable Z that is independent of a qua- 
dratic harness (Xt), taking (Xt) — {Xt + {t+l)Z) we get a quadratic harness with 
E(XsXt) = E{XsXt) + 1 + t + s + ts. To avoid such non-uniqueness, throughout 
this paper we assume 

(2.6) E{Xt) = 0, E{XtX,) = mm{t, s}. 
If (ESI) and (ESJ hold true and < s < t < u, then 

(2.7) E[Xt\T>^] = -X^ 

~ u 

and 

(2.8) E[Xt\T<,]^X,, 

see ^1 (4) and (5)]. From martingale (|2.8(l and reverse martingale H2.7(l condition 
it follows that the limits Xo \mit\^Q Xt^ and limi_^oo^t/^ exist with probability 
one. Under assumption H2.6|l wc have 

(2.9) limXi^O, limXf/i = 0, 

SO without loss of generality we may extend {Xt)t>o to include the value Xq = 
when convenient. Similarly, we may extend (^t)t>o — (^^i/t)t>o to include the 
value Xq = corresponding to i = 0. 

We now turn to the quadratic harness condition. The familiar quadratic mar- 
tingale condition associated with the martingale property (|2.8|l can be written as 

(2.10) E [X'^\ T<s] =^Xl+t-s. 
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This suggests that the quadratic harness condition associated with the simple har- 
ness property H2.3(l should be written as 

(2.11) ¥.[X^\Ts,u] =Qt,s.u{Xs,Xu), 
where 

(2.12) Qt,s.,u{x, y) = At^s.ux^ + Bt^s.uxy + Ct.s.uv'^ + Dt^s.uX + Et.s,uV + ^t,s,« 

is a quadratic form in variables x, y with time-dependent coefficients; under as- 
sumption (|2.6|l we trivially have 

(2.13) sAt^s,u + sBt^sM + uCt^s,u + Ft^s,u = t. 

Definition 2.2. A square integrable process {Xt)t>o is a quadratic harness on M_|_ 
if it satisfies conditions (|2.3|) and H2.11|l . 

The well-known examples of quadratic harnesses are the Wiener, Poisson, and 
Gamma processes. Ref. |32| identifies all quadratic harnesses with covariance 
(|2.6|) when the quadratic form on the right hand side of (|2.11l) is such that the 
corresponding conditional variance is a function of the increments X^—Xs only; this 
adds to the already listed examples two additional Levy processes: the Pascal, and 
Meixner processes. A related non-commutative form of condition (|2.11l) appears 
in [S]. The main result of [T^ asserts that quadratic harnesses with covariance 
(|2.6(l which satisfy the quadratic martingale condition (|2.10() are in fact uniquely 
determined g-Meixner Markov processes. A quadratic harness that does not satisfy 
condition (|2.1U|) is analyzed in |llj . 

2.2. Five-parameter representation. Generically, the quadratic form on the 
right hand side of H2.11|l is determined uniquely up to five numerical parameters. 

Theorem 2.2. Let (Xt) be a quadratic harness with covariance \2.(^) . Suppose that 
1^2.11]) holds with Ft^s,u ^ 0, for all < s < t < u, and that 1, Xg^ Xi, X^Xi, A"^, X"^ 
are linearly independent for all < s < t. Then there exist 77, G M, ct, r > 0, and 
q < I + 2y^aT such that 

YaT[Xt\J^s^u] = 

(2.14) (1 + 77A,,„ + eAs^u + T^l^ + tAI^ - (1 - g)A, 
for all < s < t < u, where 

(2.15) F,^^^^ iu^t)it-s) 



u(l + as) + T — qs 

is the normalizing function, and As ^u, As,u are given by 1^2. It2.^} . 

Recall that the conditional variance of X with respect to a cr-field is defined 
as Var[A:| J^] =E[X^\J']- (E[A:| J^])^ 

Remark 2.1. Time-inversion {Xt) i~> (tA"x/f) preserves the class of quadratic har- 
nesses, modifying the coefficients in (I2.11f) . More precisely, suppose (Xt) satisfies 
the assumptions of Theorem l2.2l and let Xt — tXi/t be its time inverse. Then {Xt) 

is a quadratic harness with respect to its cr-fields J^s,u ~ ^i/u.i/s, and (|2.14() holds 
with the roles of the parameters (77, 9) and {a, r) switched within each pair: 



Var Xt Ts,u /Ft 



t,S,U 
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= 1 + 0As,u{X) + V^sAX) + tAI^{X) + aAl^X) - (1 - g)A,,„(X)A,,„(X), 

where Ft s u = — . 

Ref. |18| gives criteria for time-inversion invariance of Markov processes. 

2.3. Orthogonal martingale polynomials. Suppose that a quadratic harness 
(Xt) has moments of all orders and martingale polynomials pn{x]t) of all degrees 
n > 0, that is 

(2.16) E [p„ {Xt;t)\ T<s] = pn {Xs; s) , < s < t. 

Clearly, po = 1 and pi{x;t) = x are natural initial choices, see 1)2. 8|l . Since xp„(x; t) 
is a polynomial of degree n + 1 , it follows that 

n+l 

(2.17) xpn{x;t) = y^^Ck,n{t)pk{x;t). 

k=0 

Theorem 2.3. Suppose a quadratic harness (Xt) with covariance 1^2. b]) and condi- 
tional variance \2.1^\j has finite moments of all orders and martingale polynomials 
Pn{x;t). If for each t > the random variable Xt has infinite support, then recur- 
rence i2.17\ ) holds with the infinite matrices 

Cooit) Coiit) Co2{t) Co At) ■■■ 

Cioit) Cn{t) cut) Cisit) ... 

C2l{t) C22{t) C2z{t) ... 

C^t) C33{t) ... 

dsit) ■■■ 



Ct 



given by 

(2.18) Ct=tx+y,t>0, 

and the infinite matrices X= Ci — Cq, y= Cq satisfy equation ^^\). 

From Theorem l2.3l we derive a number of equations that eventually determine the 
orthogonal martingale polynomials. Namely, it is well known that for orthogonal 
martingale polynomials Pn{x;t) recurrence 1)2.1711 holds with a tri-diagonal matrix 
C, see 115 . Writing ifTTTjl as 

(2.19) xpn{x;t) = an{t)pn+i{x;t) + bn{t)pn{x;t) + Cn{t)pn-i{x;t), n > 0, 
from ()2.18)l we get 

(2.20) Unit) = aan+it + /3„+i , 5„(t) = 7„i + Sn , c„(t) = e„t + ipn, 

and the coefficients in ()2.20)l satisfy a number of equations that result from the 
(j-commutation equation 1)1.11) . 

Setting p-i(x;t) = 0, po{x;t) = 1, pi{x;t) = x, and using ()2.6I) we see that the 
initial values are given by 

(2.21) ai = 0, /3i = 1, 70 = <5o = 0, £i = l,^i= 0. 
For n > 1, 1)1.1)1 implies that coefficients am Pmlm^m^mSn satisfy 



(2.22) 



a Ta„a„+i + (Ta„/3„+ig + (t/3„/3„+i = (Ta„+i/3„ 
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(2.23) /3 

cra,i+i(7„ + ^n+i)T + (cra„+i(5„+i + /3„+i7„)g + /?„+i((5„ + (5„+i)cr 

(2.24) /3„+ie„+l + ln5n + (yCtn^Pn 

= ((Ta„+ie„+i + + (Ta„e„)T + (cra„+i(/3„+i + 7„(5„ + 

+ (^„+l(^„+l + (5^ + /3„(y9„)o- + JnO + SnTj + 1, 

(2.25) 7„-i</?„ + (5„£„ 



(2.26) EnVn+l = £n£n+lT + En+lfuq + fnfn+lCT- 

If a„(t) 7^ for all n, then the three step recurrence (|2.19() defines a family 
of polynomials in variable x. We now show that these are indeed the martingale 
orthogonal polynomials for (Xf). 

Theorem 2.4. Suppose a quadratic harness (Xt) with covariance \2. 6|) and condi- 
tional variance \2.14^ has finite moments of all orders, parameters q, rj, 9, cr, t are 
such that (q, err) ^ (—1, 1), and the equations I^2.21}12.^) have a solution such that 
a„(i) ^ for allt>0,n> 0. Let {pn{x: <) : 7i > 0, <> 0} satisfy IFI^) and lEI^) 
for n — 0,1,... with p^i{x,t) — 0,pQ{x,t) — 1, pi{x,t) — x. Then {pn{x;t)} are 
orthogonal martingale polynomials for (Xt). 

In Sectional we give sufficient conditions in terms of parameters q,7],9,a,T for 
the assumptions of Theorem 12.41 to be satisfied, and we derive explicit three step 
recurrences for the orthogonal martingale polynomials in several cases of interest. 

Theorem 12.41 implies that if quadratic harnesses {Xt)t>o and (Yt)t>o satisfy its 
assumptions with the same parameters rj, 9, cr, t, q, then 

E{x-^xii^ . . . XI-) = ■ ■ ■ Yr:) 

for all /c > 1, < ii < t2 < • • • < tk, and tii, . . . ,nk G N. Under appropriate 
integrability conditions, this implies that {Xt) is Markov and is uniquely determined 
by the parameters 77, 9, cr, r, q. In particular, this is the case when sup„ |6n(i)| < cxd 
and sup„ |a„_i(i)c„(t)| < c», as in this case the recurrence (|2.19ll corresponds to a 
compactly supported measure, see jll Section 2]. However, the question of existence 
of such a process is non-trivial and the constructions are known only in special cases, 
see El, ini, da- 
Next, we show that the integrability assumption of Theorem l2.4l is automatically 
satisfied if err = 0. We remark that the result stated below does not use the full 
power of the quadratic harness condition, and generalizes | 1UI Theorem 2] and j32l 
Theorem 2(1°)]. Since 

Var [Xt\ J^s,u] = E [X^\ - {at,s,uXs + ?>t,.,„^«)', 
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using H2.9|l we can pass to the limit in (|2.14(l as u ^ oo or as s ^ 0. This gives 

(2.27) \&Y[Xt\T<s\ = -^-^ (dXf + r/X, + 1) , 

^ 1 + crs ^ ' 

(2.28) Var[X,|.F>J = + 0^ + 

Theorem 2.5. If a square-integrable stochastic process {Xt) with covariance \2. 61) 
satisfies ^FF^ and ^1^) with a, r > such that ar < l/24'-+io 

for some r > 2 then E{\Xt\^) < oo for all t > 0. In particular, if err = 0, then 
E{\Xt\'') < oo for all r,t>0. 

3. Proofs of the main results 
3.1. Proof of Theorem 12.21 For < x < y < z define 



^ -^y^x^z ~t~ XzBy^x,z ~\~ Z Cy^x^z U 
-^y.x,z{^ '^) 

F ' 



1 + qy,x,^ 
Vy,x,. 



F 

•^^y,x,z ~\~ zEy^x, 



We will show that the left hand sides of the above equations do not depend on 
the arguments y,x,z. To this end we prove three claims. Throughout the proofs 
of all three Claims, 0<r<s<t<u are arbitrary numbers. 

Claim 3.1. For all < x < y < z, and for f — a,T, q, rj, 9 we have 

(3.1) fy,x,z — fy^x.z^ 

provided x < y < z. 

Proof of Claim \cl.l\ Observe that from H2.3|l and (|2.11ll we get 

E[XsXt\ Tr,u] = E [X,E [Xt\Ts^u] I Tr^u] 

= at,s,«E [X^l + bt^s,uXiiE[Xs\J-rM] 
On the other hand, 

E [XsXtl J^r,u] - E [E [Xs\Tr^t] Xt\J^r^u] = a.,r,t^rE [Xt\ J^r,u] + bs,rj.^ [X^\ J'r,u] 



8 



WLODZIMIERZ BRYC, WOJCIECH MATYSIAK, AND JACEK WESOLOWSKI 



Comparing 


I the coefficients at X^, 


XrXu, Xl, Xr, Xu and 1 in 


pressions we | 


^et 




(3.2) 






(3.3) 


^t.s^u^ s,r.u 


— ^s,r.t^t,r,u ; 


(3.4) 


^t.s,u^s,r,u ~^ ^t.s.uC s^r.u 


— ^s,r,t Ct^r,u 1 


(3.5) 


^t,s,u-Ds,r,u 


— ^s,r,t-Dt,r,u 1 


(3.6) 




— ^s,r,t-^t,r,u 1 


(3.7) 




— ^s,r,t-^t,r,u • 



(In the derivation of H3.3|) we used the fact that bt^s,uO-s,r,u = o,s,r,tbt,r,u)- 
Adding and ||33I), we get 

^t,s.u (^s,r,ii ~t~ -^s,r.u ~t~ Cs,r.u 1) ~^ ^t,s,u ~t~ ^t,s,ii^s,r,iz 

Since a calculation shows that 

(3.8) 0>t^s,u ~t" bt,s,ubs,r,u ^ bs.r,t ^" ^s,r,tO>t,r,u^ 

dividing by H3.7() . we get crs,r.« = o-t,r,u which proves H3.1|l when f = a. 

We proceed similarly when / = T,q,r],9. Adding H3.2|l multiplied by r^, 
(|3.3(l multiplied by ru H3.4(l multiplied by u^, and dividing by l|3.7|l . we obtain 
Ts,r,u = Tt.r,u, after noticing that s'^at^s.u + u'^bt,s,ubs,r,u = t^bs,r,t + r'^as,r,tat.r,u- 
Equation (|3.3|) multiplied by [u—r) and divided by (|3.7I) gives qs^r.u — qt,r,u- Adding 
equations l|3.5|) and H3.6|l and dividing by (|3.7f) gives ?7s,r,ti = Vt,r,u and similarly 
multiplying (|3.5I) by r and (|3.6I) by u after dividing by (|3.7() gives 9s.r,u — Qt,r,u- 
Thus we obtained fs,r,u — ft,r,u ior f = a,T,q,ri,9. liQ<x<y<z then substitu- 
tion r = X, s = y,u = z yields (|3.1(1 for y £ (y, z); substitution r = x,t = y,u = z 
gives H3.1|l for y G {x,y), completing the proof of Claim (|3.1|l . □ 

Claim 3.2. For all < x < y < z, and f ^ a,T, q, rj, 9 we have 

(3.9) fy,x,z — fy,x,z, 

provided <x < y. 

Proof of Claim Now consider the identity 

Using f^I^ and 12.1111 wc sec that E [E [Xf\ !Fs,u] \ ^r.u] is given by 

E [At^s.uX'^ + Bt^s^uXgXu + Ct^s.uX^ + Dt^s,uXs + Et^s.uXu H- ^i,s,ti| ^ r,v\ 
= ^t,s,iiE I ^r,u\ + -St,s,u^uE ^r,ii] + Ct,s.uX^^ + -Di^s.^^E [X^l J^r,u] 

~^-Et,s,uXu ~h Et,s,u — -At^s,u{^-As.r,uXj. -\- B s_^'f^^i^Xj-X^^ -\- C g^Y-,uX^ -\- D s,r.uXf 
~^-Bs,r,uXu -h -^s,r,ti) -^t.s/u-^u ('^s.rjU-^r ^" ^s.r,uXu) Cf^s.uX^ 

~l~-^t,s,u-^s,r,u -^t^s.u- 

On the other hand, 

E [^i^l ^ r,u] — At^r,uX^ + Bt^r,uXrXu + Ci,r,ii^u "t" C'i^rjU^r + Dt^r.uXu + Ft^r,u- 
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Comparing the coefficients at X'^, XrXu, X^, Xr, Xu and 1, we obtain 

(3.10) -^t.r.u — -^t.s,u^s,r,u 

(3.11) ^t.r,u -^t.s.u^ s,r,u ~^ ^t.s,u^s,r,u 5 

(3.12) Ct_r,u — -^t-s^uCs^r-u ~(~ ^t.s,u^s,r,u H" Ct^s,u •) 

(3.13) n — -^t.s^u^ s,r,u H~ ^t.s,u^s,r,u 1 

(3.14) ^i,r,ii — -^t^s,u-Ejs,r,u H" ^t^s,u^s^r,u H" ^t^s,u : 

(3.15) Ft,r,u — -^t^s,u-^s,r,u H~ ^t.,s,u • 

Substituting the right hand sides of equation H3.10I3.15|1 for At^r,u: Bt^r,uj ■ ■ ■ , -fi,r,u 
on the right hand sides of (|3. 213 .7(1 we get 



(3.16) 








(3.17) 






— bs,r,tBt,s,u^s,r,u i 


(3.18) 








(3.19) 


D s,r.u (^^t,s,u 


~ bs,r,tAt,s,u) 


— bs,r,tDt,s,u^s,r,u ; 


(3.20) 


-^s,r,u (^t,s,u 




— ^s,r,t{J-^t,s.u^s.r.u ~^ ^t,s.,u) i 


(3.21) 


-^s,r,u (^t,s,u 


~ ^s.r^t-^t.s.u) 


— bs,r,t-^t,s.u ■ 



We can now proceed analogously to the proof of Claim mi Namely, adding H3.16|l . 
H3.17|l . H3.18|l . again taking into account 1)3.8(1 . and dividing by ((3.21(1 . we get as,r,u = 
(Jt.s,u- Adding ((3.1()() multiplied by r^, 1(3.17(1 multiplied by ru and 1(3.18(1 multiplied 
by , and dividing by ((3.21() . we obtain Tq j-^h — '^'t s u- 

Equation 1(3.17(1 multiplied 

by {u — r) and divided by 1(3.21(1 gives qs,r,u — <lt,s,u- Adding equations ((3.191) and 
((3.20(1 and dividing by 1(3.21(1 gives rjg^r.u — Vt,s,u and similarly multiplying ((3.19(1 
by r and 1(3.20(1 by u after dividing by ((3.21(1 gives 9s.r,u = Ot^s.u- Thus we obtained 
fs,r.u = ft,s,u for / = a,T,q,r],d. By Claim EH] /s,r,u = ft,r.u, so Jt.r,u = ft.s,u, 
ending the proof of Claim □ 

Claim 3.3. For all < x < y < z, and / = cr, r, q, rj, 6 we have 

(3.22) fy,x,z — fy,x,z: 

provided z > y. 

Proof of Claim UTM This follows from Claim bv the time- inversion {Xt) i-^ 
(tXi/t). Alternatively, one can repeat the previous arguments, starting with the 
identity 

E Tr;r,] =E[E[X^\ Tr^t] \ ^r,n] ■ 

□ 

Conclusion of proof of Theorem \2.'A Now, it is easy to deduce that ((3.1(1 . 1(3.9(1 and 
((3.22(1 imply that functions a, r, q, rj, 6 are in fact constants. Indeed, given < xi < 
yi < z\ and ^ < xi < yi < zi, and / = a,T,q,r],9, with x :— min{a;i,a;2} and 
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z := maxjzi, 2:2}, we see that 

fyuxi.zi = fyuxuz by ClaimESl 

= by Claim O 

= fy2,x,z by ClaimEm 

= fy2,x2,z by Claim EiS 

= fy2,x2,z2 by ClaimEiS 



(If X2 = X or Z2 — z, one or both of the last two steps is unnecessary.) 
So, for all < s < t < u we have 



(d.2jj (7 = — — — — , r 



Ft Ft 
(3.24) l + q= ,7?- -,9^ . 

^t^s.u ^t^s.u ^t^s.u 

The above equations, along with H2.13|l . make a system of linear equations in vari- 
ables At^s.u, ■ ■ ■ , Ft.s,u- This system must be solvable for all < s < t < u, so its 
determinant {u — s)"^(u(l + as) + r — qs) 7^ 0. As w — > 00, s — s- 0, the expression 
u{\ + as) + T — qs is positive. Thus u(l + as) + t — qs > 0. Taking the limits 
s — > 0, u — !■ 0, s = 1, u — !■ cx), and s —> u, we get r > 0, a > and g < 1 + au + t/u 
respectively. Minimizing the latter over u > we get q < 1 + 2^/aT. 

The unique solution of the system of equations is given by H2.15|l , (|3.25I3.29|I . A 
calculation shows that l|2.14(l holds. □ 

Remark 3.1. A calculation shows that formula H2.14|l is equivalent to the following 
formulas, valid for all < s < i < w 



(3.25) At,s,u 

(3.26) Bt.,,u 

(3.27) Ct^s,u 

(3.28) A,s,« 



(3.29) E, 



't,S,U 



{u 


-t)[u{l + 


at) + 


T-qt] 


{u 


-s)[u{l + 


as) + 


T-qs] 




(u - t)[t - 


-s){l- 




{u 


-s)[u{l + 


as) + 


T-qs] 


it 


-s)[t{l + 


as) + 


T - qs] 


{u 


-s)[u{l + 


as) + 


T-qs] 




{u-t){t- 


s){uri 


-0) 


{u 


-s)[u{l + 


as) + 


T-qs] 




[u~t){t- 


s){e- 


' sv) 


{u 


-s)[u(l + 


as) + 


T - qs] 



and Ft,s,u is given by H2.15|l . 

3.2. Proof of Theorem lOl For non-commutative variables X,y, let 
(3.30) Qt,,,„(x, y) = At,s.u>^^ + Bt.s,u>^V + Ct,s.uf + A,s,«x -f £:t,,,„y -|- Ft 



t.s^u , 



recall the quadratic form H2.12|l in commuting variables x,y. (Note that the order 
of X, y at Bt.s,u is changed when compared to its dual version (|4.35|) .) 
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Lemma 3.4. Under the assumptions of Theorem \2.'A 

(3.31) Ct — at^s.uCs + bt^s,uCu , 

(3.32) — Qt^s,uiCs,Cu) , 

where the coefficients in fy.c/U\) are given by f2.15\) and j^*. 25t!^. 2y\) . 

Proof. For a polynomial iy9 : R — > M consider vectors 

v,,t = [Ei^{XMXt-,t)),E{ipiXMXt-,t)), . . . ,E{ipiX,)pkiXt-,t)), . . .] eR^ . 

Since for s > random variable Xs has infinite support, polynomials 1, Xg, X^, . . . 
are linearly independent, and the corresponding orthogonal polynomials are non- 
degenerate. This implies that as we change (p for a fixed s > 0, vectors of the form 
Vg^s are dense in R°° , equipped with the product topology. Indeed, applying the 
Gram-Schmidt orthogonalization process to Po{x; s),pi{x; s), . . . we get a sequence 
of orthogonal polynomials qo,qi,... such that M{qi.{Xs)-^) ^ 0. Therefore, for 
any w = [wq,Wi, . . .] G M°°, and n > we can find numbers uq, Wi, . . . , u„ such 
that the first n coordinates of v^^s corresponding to <y9 = X]j=o ^i'^i ^'"^ equal to 
Wo, wi, . . . ,w„. 

This implies that to verify the identities, it suffices to verify that the identities 
hold true when multiplied from the left by Vs,s. 

By the martingale polynomial property, E{(p{Xs)pk{Xs; s)) — E{Lp{Xs)pk{Xt;t)). 
Therefore, from H2.17|l it follows that Vg.s = 'Vs,t and 

Vs^s xCt= Vs,t X Cf 

E ^^{x,)Y,c,At)pAXt;t)^ ,E ^^{x,)Y,c,At)pAXut) 

= [Eiip{X,)XtPo{Xt;t)), E{ip{X,)XtPi{Xt;t)), . . . ,E{AX,)XtPk{Xt;t)), . . .] 
= [E{ipiX,)XtPoiXu; u)),Ei^iX,)XtPi{Xu; u)),..., E{ip{X,)XtPk{Xu; u)),...] 
= at,s,u [E{ipiX,)X,pk{Xu;u)) : fc > 0] + [E(^(X,)X„pfc(X„; u)) : fc > 0] . 
Using the martingale polynomial property again, we see that v^.s x Ct equals to 
at,,,„ [E{ipiX,)X,pkiX,; s)) : k > 0] + bt,s,u [E(^(X,)X„pfc(X„; u)):k> 0] 

0't,s,u^s,s X Cg + bt.s^u^ s,u X C,( — Wg^g X (ai^s^uCs + &f,s,uCu), 

proving (|3.31() . Similar reasoning proves (|3.32() : 

V,,, X = v,,t X = [E{p{Xg)XtPk{Xut)) : fc > 0] X Ct 
= [E{p{Xg)X-ipk{Xt-t)) : fc > 0] = [E{p{Xg)X-^pk{X^-u)) : fc > O] 
= ^{ip{Xg)E Pfc(X„; u)):k> O] 

= [E{^{Xg)Qt^g^AXg,Xt)pk{X^-u)) : fc > 0] 
= At^g.u \E{ip{Xg)Xlpk{Xg-s)) : fc > 0] -I- Bt,g.u [E{ip{Xg)XgXuPk{Xu-,u)) : fc > 0] 
+Ct,g,u [E{^{Xg)Xlpk(Xu-u)) : fc > 0] + Dt,g,u [E{ip{Xg)XgPk(Xg- s)) : fc > 0] 
+Et,g,u [E{^{Xg)XuPk{Xu]u)) : k > 0] +Ft„,,„ [E(</j(X,)pfe(X,; s)) : fc > 0] 
= At,g.u^g^g X C2 + Bt.g^^ [E{^{Xg)XgPk(Xg; s)) : fc > 0] x C„ 
+C't,s,iiVs,s X -|- Dt^g,u'Vg^g X -f- Et g ^Vg g X C„ + Ft g y^-Vg g 

+ Ct. s,uC„ -|- Dt.s^uCg + Et.s^uCu + Ft,s^ui) ■ 

□ 
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Lemma 3.5. Suppose a linear operator C{t) = tx + y satisfies relation IS.H'^) for 
some quadratic form Qt,s,u given by \S.S(]\) . If the coefficients of Qt.s,u are given 
by fin. 15}) . then Hl.l}) holds true. 

Proof. Expanding the expressions like = t^X^ + t{xy + yx) + and using (|3.32() 
we get 

(3.33) t^x^ + t(xy + yx) +f = At,sAs'^>^'^ + ^(xy + yx) + y^) 
+Bt,s,u[f + + sxy + suy?) + Ct^^,„(u^x2 + u(xy + yx) + y^) 

HDt,s,u + Et,s,u)y + {sDt,s.u + uEt^s,u)y- + Ft,s,J- 

Applying the relations p. 13(1 . H3.23|l . and H3.24|l to the coefRcients at each of the 
monomials in (|3.33() we get 

xy + yx = TX^ + cry^ + r;y + 6'x + I + (1 + q)yx, 

which is the same as □ 

3.3. Proof of Theorem 12. 4L First we shall show the martingale property H2.16|l . 
To this end, we proceed by induction on n. Trivially, H2.16|l holds for n = 0, 1, see 
(|2.8|) . Assume then that n> 1 and H2.16|) holds for 0, 1, . . . , n - 1 and aU < s < i. 
We start from calculating E [XtPn-i (Ar„; u)\ J-<s] in two ways. On one hand 

E [XtPn-l {Xu, u)\ .F<,] - E [E [Xt\:F,^u]Pn-i (Xu, u)\J^<s] = 
= E[{at,s,uXs + bt.s,uXu)Pn-i iXu;u)\T<s] = 

= at,s,uXsPn-l (Xs] s) + bt^s.vM [XuPn-l {Xu, u)\ J^<s] = 

= at,s,u{an-lis)pn (Xs; s) + bn-l{s)pn-l (Xs] s) + Cn-lPn-2 (Xs', s) ) + 

+ bt,s.u{a„-i{u)E [pn (X„;m)| T<s] + 6„_i(m)p„_i (X^is) + c„_i(m)p„_2 (Xs^s)). 
On the other hand, 

E [Xtpn-i (A„; u)\ J^<s] = E [XtE [p^-i (X„; u)\ J^<t]\ J^<s] = E [XtPn^i {Xt;t)\ T<s 

= a„_i(i)E [pn {Xt\t)\T<s\ + bn-i{t)pn^i {Xs;s) + c„_i(i)p„_2 (A^; s) . 

Thus comparing the right hand sides of the above equations we obtain the fol- 
lowing equation: 

bt^s.uan-iiu)^. [pn (A„; u)\T<s] - a„_i(t)E [p„ (Af ; t)| T<s] = 
= -pn [Xs; s) at^s,uan-i{s) +Pn-i (As; s) {b.n-i{t) - at^s,ubn-i{s) - bt,s,ub„-iiu)j 

+Pn-2 {Xs] s) (^C„_i(t) - at^s,uCn-l{s) - &t,s,uC„-l (m)^ . 

A trivial verification, using H2.2U|I . shows that 

bn^i{t) ~ at^s.uK-lis) ~ bt^s.ubn-l{u) = 0, 
Cn-l{t) - at,s,uCn-l{s) - bt,s,uCn-liu) = 0. 

Hence denoting 

A- E [p„ (A„; u)\ J^<s] , 3^ := E [p„ {Xt;t)\ T<s] , 

we have 

(3.34) bt,s,uan-iiu)X ~ an-i{t)y = -at,s,ua„-i{s)pn (A^; s) . 
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To obtain a second equation for X and y let us consider E [Xfpn-2 {Xu, u)\ J^<s] ■ 
On one hand 

E [Xfpn-2 (^«; u)\ ^<,] = E [XfE [p„_2 u)\ J^<t]\ J^<s] - 

= E[xfpr,-2 {Xt;t)\T<,] . 
Setting ak{t) — bk{t) — Ck{t) = for fc < 0, a repeated application of H2.19|l gives 
(3.35) 

x'^Pn^2 ix;t) = pn ( x; t) a„_2 (<)a„- 1 (i) + Pn-i {x;t) a„_2(<) [bn-2{t) + &„-i(i)] + 

+ Pn-2 ix;t) [a„_2(^)c„-l(0 + bf^-2it) + Cn-2it)an-3{t)] + 
+ Pn-3{x;t) Cn-2{t) [bn-3.{t) +fc„-2(i)] +^,1-4(2;;*) C„-3(i)c„-2 

we obtain 

(3.36) E [X'^Pn-2 {Xu;u)\ T<s] = a„-2(t)an-i(i) 3^+ 

+ Pn-1 [Xs] s) a„_2(^) [K-2{t) + hn-l{t)] + 
+ Vn-2 {Xs] s) [an-2{t)Cn-l{t) + 6^-2(i) + C„-2 (i)ari-3 (0] + 
+ p„_3 {Xs-, S) C„_2(i) [hn-z{t) + hn-2{t)] + 

+ Pn-4 {Xs; s) Cn-3{t)Cn-2{t)- 

On the other hand, one can rewrite E [Xfpn-2 {Xu, u) \ ^<s\ as 

E [Xfpr,-2 {Xu, U)\j^<s] = E [E [X^'l Pn-2 {Xu: u)\ .F<,] 

= At,s.uX^Pn-2 {Xs-, S) + Bt,s,uXsE [X„p„_2 {Xu]u)\ T<s\ 
+Ct,,,„E [Xlpr,-2 {Xu, u)\T<s] + Dt,s^uXsPn-2 {Xs] s) 
+Et,s.vM [XuPn-2 {Xu, u)\T<s] + Ft^s.,uPn-2 {Xs] s) . 

After some algebra one gets 

(3.37) E [Xfpn^2 {Xu;u)\ T<s] = X Ct,s,ua„-2{u)an-l{u) 

+ Pn {Xs; s) an-l{s)Rn{t, S,u) +Pn-1 {Xs; s) Sn{t,S,u) + Pn-2 {Xg; s)Tn{t,S,u) + 

+ Pn-3 {Xs; s) Un{t, S, u) + Pn-4 {Xs; s) Vn{t, S, u) , 

where 

Rn{t,S,u) = At^s,uan-2{s) + Bt^s,uan-2{u), 

Sn{t, S, u) = At^s,uan-2{s) [&«-2(s) + fcn-l(s)] + Bt^s.u [a„_2 (")&«- 1 («) + a„_2 (s)6„_2 (")] 
+ Ct,s.uan-2{u) [bn-2{u) + 6„-l(u)] + A,s,«&n-2(s) + Et^s .ubn-2{u) + Ft,s,u , 

T„{t,S,u) = At,s.u[a„-2{s)Cn-l{s) +bl_2{s) + Cn^2{s)an-3{s)] 

+ Bt,s.u [a„_2(w)c„-l(s) + bn-2{u)bn-2{s) + C„-2 (M)a«-3 (s)] 
+ Ct,s.u [a„_2(")c„-l(M) + bl_2{u) + Cn-2{u)an-3{u)] 
+ Dt^s,ubn-2{s) + Et^s.,ubn-2{u) + Ft.s,u , 

Un{t, S, u) = At^s,uCn-2{s) [bn-3{s) + bn-2{s)] + Bt.s,u [&n-2 (u)c„-2 («) + C„_2 ('«)fcn-3 (s)] 
+ Ct.s,uCn-2{u) [fo„-3(") + bn-2{u)\ + Dt^s.,uCn~2{s) + -Bt^s,«C„_2 (u) , 

Vn{t, S, U) = At^s,uCn-2{s)Cn-3{s) + Bt,s,«C„_2 (w)c„_3 (s) + Ct,s,„Cn_2 (w)c„_3 (m) . 
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Comparing the right hand side of H3.37|l with the right hand side of H3.36|l we 
get the second equation with unknowns X and y: 

(3.38) Ct,sMan-2{u)an-i{u) X - an-2{t)an-i{t) y = -pn {Xs] s) Rn{t, s,u) 

+ Pn-1 {Xs; s) (^a„_2(i) [bn-2{t) + bn-l{t)] - Sn{t, S, u)^ 
+ Pn-2 [Xs] s) {an-2{t)Cn-l{t) + &^_2(^) + C„_2 (OOn-S (0 " Tn{t,S,u)^ 
+ Pn-3{Xs;s) (^Cn-2{t) [bn-sit) +fori-2(^)] - Un{t,S,u)j 

+ Pn-i {Xs;s) (^C„_2(i)c„-3(0 - Vn{t,S,u)y 

A calculation based on (|2.2Q() and l|3. 2513. 29(1 gives 

a„_2(t) [bn-2it) +bn-l{t)] ~ Sn{t,S,u) 

it-s){u-t) 



I'T'n — 1 ~r (7T(y.n — l wn — 2 ~r '~fn—l ) 

u{l + as) + r — L 

+ 0-/3„-l {6,1-2 + Sn-l) + q iPn 



an-2(i)c„-l(t) + b'^n-2{t) + Cn-2 (i)an-3 (0 " Tn{t, S,u) 



u[l + as) + T — qs I 

+ g'(7n-2^n-2 + ) + CT ((5^_2 + Pn-2(l)n-2 + /3„_10„-l) 



c„-2(i) [bn-ait) + bn-2{t)] - i7„(i, S, u) 



(i - s)(u - t) 



u{l + (Ts) + r — gs 

+ Te„_2 (7n-3 + ln-2) - &n-2£rL-2 + V4'n-2 - 0n-27n-3 
+ a {Sn-3 + Sn-2) 0n-2 + Q {S 



c«-2(t)c„-3(i) - ^(i, 

(i - s)(w - t) 



-3 (Te„-2 - 4'n-2) + 4>n-3 {q£n-2 + 0-(t>n-2)] , 



u{l + as) + T — qs 

Therefore, by H2.23|) - H2.26|) the coefficients at Pn-i {Xs;s) , . . . ,Pn-i [Xg] s) on the 
right hand side of (|3.38|l vanish, and (|3.38() is equivalent to 

(3.39) Ct.s,uan-2{u)an-i{u) X - a„_2(<)an-i (0 y 

= -an-i{s)pn {Xs^s) (At^s,„a„_2(s) + Bt,s,«a„_2 («)) • 
Subtracting from (|3.39|) equation H3.34|l multiplied by a„_2(^), we get 

(3.40) a„_i(u) {Ct^s,uan-2{u) - fot,s,„a„-2(t)) X 

= -ari-l(s)Pn [Xs] s) {At,s,uan-2{s) + Bt.s ,ua,n-2{u) - at,s,uan^2{t)) ■ 
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Now we notice that 

Ct,s,uO'n-2{u) — bt^s,uO"n-2{t) = 

{t-s){u-t) 



[u — s) [u{l + as) + T — gs] 
see ()4.2|l and (|2.15|l . Similarly, 

{t-s){u-t) 



[{qs ~ T)(jan-i + (1 + so-)/3„_i] 



(m — s) + as) + T — qs] 
Since (|2.22|) implies that 



[{qu - T)aan-i + (1 + ua)l3n-i 



(3.41) a„_i(u) ((gs - T)aan-i + (1 + so-)/3„-i) 

= a„_i(s) ((gu - r)(Ta„_i + (1 + ua)Pn-i) , 

and a„_i(u) 7^ by assumption, therefore, (|3.4U|) becomes 
(3.42) 

{{qs - T)aari-i + (1 + S(t)/5„_i) X = ((gs - r)crQ!„_i + (1 + scr)/?„_i) p„ (X^; s) . 
It remains to verify that 

(3.43) {{qs - T)aa„^i + (1 + sa)/3„_i) ^ 0. 

Since (3n-i — a„_i(0) 7^ by assumption, this is trivially true if a = 0. Suppose 
a > and l|3.43ll isn't true. Since a„_i(u) 7^ from H3.41|l we see that then the left 
hand side of (|3.43() must vanish also when s is replaced by u. Thus 

(3.44) a{an-iq + /3„-i) = 0, /3„_i - aran-i = 0. 

Since /3„_i 7^ 0, the determinant of the system of equations (|3.44|) is zero, q+ar = 0, 
and the second equation gives ar ^ 0. Since oi\ = 0, from (|2.22() with q = —ar 7^ 
we get aj = (ij for aU j > 1. As err 7^ 1, this contradicts 13.44|l . Thus H3.43(l holds. 
By (E33Il, from we get 

X = pn {Xs;s), 

which means that {pn{x] t)}oQ is a martingale for every n > 0. 

Now we shall show that polynomials {pn{x; t)}„ are orthogonal with respect to 
the distribution of Xt. Since the polynomials satisfy the three term recurrence 
(|2.19|l . it suffices to show that for all fc = 1, 2, . . . 

(3.45) E{pk{Xt;t))^0. 

Clearly, by the martingale property, /i^ = K{pk {Xt; t)) does not depend on t. 

We proceed again by induction. From 14.2|l we see that ai{t) = 1 + at so 
(1 + at)p2{x; t) = — bi{t)pi{x; t) — t. Since pi{x; t) = x, assumption H2.6|l implies 
that holds for fc = 1, 2. 

Assume then that for some n > 3 we have /ifc = for all 1 < k < n — 1. By 

E{XtPn {Xt;t)) = a„(t)E(p„+i {Xt;t)) = a„(i)//„+i. 
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By ()2.7(l . the martingale property of (p„(x;i))„, (|2.19(l and the induction assump- 
tion, the left hand side of the above equation transforms into 

EiXtPn{Xt;t))^E(^^E[X,\T>t]pniXut)') = ^e(x,E [p„ (X*; t)| .F<,] ) = 

= -EfXsPn {Xs; s) ) = -a„(s)E(p„+i (X^; s)) = -a„(s)/i„+i, 
s \ / s s 

so 

'a„{s) Unit) 



t 



Pri+1 = 0. 



< ( Cilx" + C2/X + ) (Pr(|X| >x)+ Pr(|r| > x)) . 



Since 

s t ts 

we get /i„+i = E{pn+i {Xt]t)) = 0, and the orthogonality of {Pn(a;;i)}„ is proved. 

3.4. Proof of Theorem EISJ 

Lemma 3.6. Suppose X,Y are square-integrable standardized, and there are con- 
stants e, j4, _B > 0, < p < 1 such that 

(3.46) E[{X ~ pYf\Y] < A + B\Y\ + ^^^^j-^^Y'^, 

(3.47) E[{Y - pXf\x] < A + B\X\+ ^^~f^f X^- 

If 2e < p < 1, then there are constants Ci, C2 < cxd such that for x > we have 

(3.48) Pr(|X| > 2a;/p) +Pr(|y| > 2a;/p) 

4£2 1 + 

p2 + £2 1 - 

Proof. Enlarging e if necessary, without loss of generality we may assume that 
£ > 0. 

Let N{x) = P{\X\ >x)+ P{\Y\ >x),K = 2/ p. Throughout the proof, Ci, C2 
denote positive constants which might differ at each occurrence. 

The event {\X\ > Kx}, where a; > is fixed, can be decomposed into the sum 
of two disjoint events {\X\ > Kx} n > x} and {|X| > Kx} n < x}. 
Therefore, denoting 

Pi (a:) = Pr(|X| >a:,|y| >a:), 
P2{x) = Pr{\X\>Kx,\Y\<x), 
Psix) - Pr{\Y\>Kx,\X\<x), 
we get Pr(|X| > Kx) < Pi{x) + P2{x), and hence by symmetry of the assumptions 

(3.49) N{Kx) < 2Pi(x) + P2{x) + P^ix). 
To estimate Pi (a;), we observe the following. 

Claim 3.7. If {\X\ > x, |y| > x\, and < a < 1, then either 
{X - pYf > a(l - pfY^ + (1 - pf{l - a)x^ 

or 

{Y - pXf > a(l - pfX^ + (1 - pf{l - a)x\ 
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Indeed, suppose that both inequalities fail, i.e. on the set {\X\ > x, \Y\ > x} we 
have {X - pYf < a(l - p)^X^ + (1 - p)^{l - a)x^ and (Y - pX)^ < a{l - pfY^ + 
(1 — /o)^(l — a)x'^ . Adding the inequalities we obtain 

(1 + p^){X'^ + y^) - ApXY < a(l - pf{X^ + + 2(1 - pf{l - a)x2. 

Since ipXY < 2p{X^ +Y^), this gives 

(1 - - a){X^ + Y^) < 2(1 - p)\l - a)x\ 

However, since \X\ > x, \Y\ > x, we have 

(1 - pfil - a){X' + Y^) > 2x2(1 - p)\l - a), 

a contradiction. 

Claim im with a = 1/2 implies 



Pi{x) < Pt{\X - pY\ > (1 - p)VY^ + x^/V2, \Y\ > x) 



+ Pr{\Y - pX\>{l^ p) VX2+^/V2, \X\ > x). 
From conditional Chebyshev's inequality and l|3.46|l we get 



Pr{\X - pY\ > (1 - p)\/Y^ +x^/V2, \Y\ > x) 
< 2 / A + B\Y\ + il-p')e^Yyil + e') ^^ 



\Y\>i 



(1 -p)2(r2 +2-2) 



l\Y\>x P)^X^ J\Y\>x - P)^\Y\x J\Y\>x {l-pfY'^ 

< Pr(|r| >x) + — Pr(|r| >x) + :^—!-Pi{\Y\ > x), 

where Ci = 2A/(1 — pY , C2 ~ 25/(1 — p)^. Since the assumptions are symmetric 
in X, y, this shows that there are constants Ci, C2 < 00 such that 



a; 1 + £^ 1 — p 

To estimate P2{x), we use the trivial estimate \Y — pX\ > p\X\ — \Y\, which shows 
that the event {|X| > Kx, \ Y\ < x} implies 



\Y - pX\ > e\X\ + {{p -e)K~l)x^ e\X\ + ^ — —x 



Therefore 



P2{x) < Pr(|r - pX\ > e\X\ + ^ — —x, \X\ > Kx) 

P 

= [ Pr(|y - pX\ > e\X\ + ll^^x\X) dP 

J\X\>Kx P 

< f A + B\X\ + {l-p')s'xyil+e') ^^ 



\X\>Kx 



ie\X\ + ip-2e)x/pr 



\X\>Kx 

{{p-2e)x/pY 

J\X\>Kx 

e\X\{p-2e)x/p 

\x\>Kx{l + e^)e^X^ ■ 
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Since N{Kx) < N{x), this shows that 

X 1 + £^ 

where Ci = p^A/ (p — 2e)^, C2 = pB / {e{p—2e)). By symmetry of the assumptions, 
we also have 

X'^ X 1 + £^ 

Combining H3.49|l . H3.50|l . H3.51|l . and (|3.52|) we obtain that there are constants 
C^^Ci > such that 

N{Kx) < ; ' + ^ ' + ^N{Kx) + ^- -N{x), 

which imphes H3.48|l . □ 



Lemma 3.8. Suppose X,Y are square-integrable standardized, 'E{\X\p)+E,{\Y\p) < 
00 for some p > 1 and the assumptions of Lemma \S.f)\ are satisfied with constants 
£>0,0<p<l such that 

(3.53) 2P+3£2i±^ < pP+\ 

1 - P 

Then E(|X|P+i) + E(|r|P+i) < 00. 

Proof. Clearly, (|3.53|) implies 2e < p. Indeed, 

1- p 

We use Lemma [3.()l and we use the notation N{x) — P{\X\ > x) + > x), 

K = 2/ p introduced in its proof. Fix Af > 0. Then, noticing that K > \ from 
Lemma [3.61 we get 

{p+l)xPN{x)dx = KP+'^ / {p+l)xPN{Kx)dx 



<^—Ci {p-l)xP-^N{x)dx + ^—C2 pxP-^N{x)dx 
P - 1 Jo P Jo 

+ P^P+i / (^p+i)j;PN{x)dx 



p^ + 1 — p 

< P±lci{Ei\X\P-') + E(|rri)) + P±lc2mx\P) + Ei\Y\P)) 
p-1 p 

+£^2P+^ / pP+'^ / {p+l)xPN{x)dx. 
pJa 

Therefore, if (|3.53l) holds true, then supj^/^p jQ \p+i)xPN{x)dx < 00 which implies 
E(|X|P+i) +E(|r|P+i) = J^{p+l)xPN{x)dx < 00. □ 

We will also need the following simple observation. 

Lemma 3.9. Fix r > 1. Under the assumptions of Theorem \2.5\ if E{\Xt\'^) < cxd 
for some t > 0, then E(|Xt|'') < 00 for all t > 0. 

Proof If E{\Xt\'') < 00 for some t > then E{\Xs\'') < 00 for all s < i by lfTH|l . 
Similarly, E(|X„|'') < 00 for all u > t by (ET)). □ 



QUADRATIC HARNESSES 



19 



Conclusion of vroof of Theorem \2.5i Let r > 2 be fixed. Since our proof relies 
only on inequalities, if ctt = we can increase slightly the values of parameters 
cr, T on the right hand sides of (|2.28|) and (|2.27|) to ensure <7',t' > and a'r' < 
2-4r-io^ replacing both equalities by the appropriate inequalities; thus without loss 
of generality we may assume that ar > 0. 

Let p = r — 1 > 1. Define — ^far jl, to — t/e^^Sq = e^/a. Then the 
correlation coefficient between Xs,, and Xt^ in p — \/ so/to — e'^/^/ar = 1/2. 
From the inequality form of 12.28|l and (|2.27l) we deduce that random variables 
X = XsJ^ and Y = XtJVto satisfy ^^1^ and (|CT7| with 

A = i-4, B = 1-4 max{h|e/V^, \0\V^/s}. 
1 + 1 + 

Condition (|3.53|) is satisfied since with r = p + 1 we have (arY^^ < 1/2'^p+'^. Thus 
1 - P 

Therefore, by Lemma f3. 81 IE(|Xt(,|'") — E(|Xt„ 1^*+^) < oo and therefore by Lemma 
El E(|Xt|'') < c» for alH > 0. □ 



4. Examples of explicit recurrences 



In this section we will say that (Xt) is a quadratic harness with parameters 
q,r],9,a,T if H2.3|l . H2.6(l . and l|2.14|l hold. Our goal is to derive explicit versions of 
()2.19|l at the expense of additional assumptions on the parameters. 

Theorem 4.1. Suppose (Xt) is a quadratic harness with parameters < err < 1, 
— 1 < g < 1 ^ 2y^CTT. Moreover, assume that for each t > Q random variable Xt has 
moments of all orders and infinite support. Then (Xt) has orthogonal martingale 
polynomials which are given by recurrence 1^2. with po(x,t) = 1, pi{x,t) — x, 
where the coefficients of the recurrence are linear functions of t, 

(4.1) an{t) = aan+it + /3„+i , 6„(f) = 7„< + (5„ , c„(t) = (/3„i + ra„)a;„ , 

which are determined as follows. 
(i) The initial conditions are 



ai = 0, /3i = 1, 70 
(ii) Sequences (q!„), (/?„) satisfy 



0, oji = 1. 



(4.2) 



Otn+l 




' q 1 ■ 


X 


an 






—err 1 





n> 1. 



Moreover, denoting 

An,fc = Pnf^n+k — (^Ta-aOLn+k 

we have Xn,k > for all n> 1, k > 0. 
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(iii) Setting Ao,2 = ao = Po = 0, sequences (7n), ((5n) satisfy the following 
system of linear recurrences 

7r! + l = T— (A„.27n + {^,1+2^,1 " Pn+2an)crSn) 

(4.3) +^^(^™„+, + 0(3^^,) + p±^^eaa^+^ + 77/?„+i), 

'^n+2,0 

<5n+l = Y~ (An,2'5„ + (q;„+2/3„ - (3n+2a.n)T^n) 

^ri+2,0 

+ /^"^^ (??TQn+i + 6'/3„+i) + ™"+^ (6>ga„+i + ?//3„+i), n > 0. 

-^n+2,0 '\n+2,0 

(iv) Sequence {ujn) satisfies the linear recurrence 

(4.4) UJn+l = 7 

, 1 + 7«(T7n - +0) + 5n{qin + oSn + f]) ^ o 

H T , n>2 

with the initial term 

'"^ - + ''^ {1 - arni - ar{2 + q)) ' 

Remark 4.1. If —1 < 9 < 1 — S-^ctt, then the exphcit solution of recurrence 
is 



(4.5) a„ = ^ r— , A 



A+ - A_ ' ' " A+ - A_ 

where 



A± = - (l + (7±v/(l + 9)'-4(g + aT) 
The proof of Theorem 14. II rehes on Theorem 12.41 and the foUowing lemma. 
Lemma 4.2. If ar < 1 and — 1 < ? < 1 — 2y/aT, then the solutions of i4-^ satisfy 

Pn > \foTan > 0. 

In particular, 

PnPn+i > <7Tanan+i, (3^ > (TTaf^ , for all n>l. 

Proof. From the explicit solution in Remark |4. II we see that |A_ | < |A+| and hence 
a„ > for n > 1. We now prove by induction that P„ > ^foToin for all n. This is 
trivially true for n = 1, see (|2.21(l . Suppose the inequality is satisfied for some n > 1. 
Then (|4.2|l implies that /3n+i — ^faroin+i = (1 — \faT){in ~ \/oT{q + ^/GT)an > 
(1- V^)(/3„- V^a„) >0. □ 

Proof of Theorem \4-l\ By Lemma [4.21 a„(<) = aan+it + (3n+i > Pn+i > for all 
t > 0, n > 0. Moreover, recurrences H4.,S|I . H4.4(l are well defined and have unique 
solution. Therefore, to end the proof we need only to verify that with 

(4.6) e„ = UJ„l3„, (pn = TUJuO-n, 

the equations H2.21I2.26|I are satisfied. Trivially, 12.21(1 holds true. It is easy to 
verify that the solution of (|4.2(l satisfies H2.22|l . Using (|4.6|l from (|4.2(l we also get 
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It remains to show that equations (|2.23(l , (|2.24(l , and p. 25(1 hold. We will rewrite 
these equations using 1(4.2(1 and its equivalent form 

(4.7) qf3n+i + crran+i = [q + crT)f3n, Pn+i - an+1 = -(q + (TT)an, n > 1. 
Equation ((2.23(1 is equivalent to 

7n+i(/3n+i - crran+i) - Sn+ia{(3n+i + an+iq) 

= ((TTa„+i + /5„+ig)7„ + (5„(t(/3„+i - a„+i) + aa„+id + Pn+iV- 
Using 1(4.2(1 and 1(4.7(1 . we can rewrite it as 

(4.8) 7„+i/3„+2 - o-5„+ia„+2 

= {q + (7T)finln - (? + (7T)aan5n + CFUn+lO + Pn+lV- 

Similarly, ((2.25() for our sequences rewrites as 

(4.9) UJn {5n+l{Pn+l " CrTa„+i) - T^fn+liPn+l + qa,i+l)) 

= UJn {{qPn+i + <JTan+i)dn + r(/3„+i - a„+i)7„ + 0f3n+i + rjTan+l) ■ 

Both of these equations are satisfied; in fact, ((4.3() was obtained by solving the 
system of equations ((4.814.9(1 when > 0. 

To verify that 1(2.24(1 holds, we substitute ((4.6(1 : using ((4.2(1 and Lemma we 
get 

(in{qPn + tTTa„) + (JTaniPn ~ ^n) 

^n+1 — 7, Q 

Pn+lPn+2 - CrTa„+l 0,1+2 

_^ 1 + InjT^n -5n+0)+ Snjqjn + CT'^k + V) 

Pn+lPn+2 - CrTan+lOln+2 

Since ((4.7(1 holds, thus ((2.24() is satisfied, too. □ 
Remark 4.2. From ((4.3(1 we get 

(4.10) ^,^!^,^,-^- + ^ 



1 — (TT ' 1 — err 

which implies that the initial recurrences are 

(4.11) xpo{x]t) ^ pi{x;t), 

(4.12) xpi{x-t) = (1 + (Jt)p2{x; t) + ^ ^""f ^ ^"^^ ^ ^\ i{x] t) + tpoix- 1). 

1 — err 

Thus the first three orthogonal martingale polynomials are 

(4.13) pQ{x]t) ^l, pi{x]t) ^ X, 

1 2 {rj + ea)t + r^T + e ^ t 



P2{x;t) = -X — —X 



1 + CTt (1 - crr)(l + crt) 1 + crt 
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4.1. Free quadratic harnesses. Free harnesses have parameter 

(4.14) q = -o-r. 

The adjective "free" is motivated by the fact that when ar = this choice is related 
to free convolutions that arise in free probabihty, see |12[ Section 4.3] and |lll 
Section 4]. In general, from Theorem l4.1l it is clear that this choice of q significantly 
simphfies the recurrences. It is not obvious whether this case is further related to 
free convolutions, see however |22| . 

Assuming ar < 1, H4.14|l imphes —l<q<l — l^far^ so it foUows from Theorem 
14.11 that the orthogonal martingale polynomials exist. It is easy to check that the 
solution of (|4.2|l is 

(4.15) a„ = (1 - fTT)"-2, /?„ = (1 - ar)"-2^ n > 2. 
For our choice of q and ri > we have 



Using this identity, from H4.3|) we get for n > 2 

77 + 20CT + ?;(TT 

(1 - or) 



^ 6* + 2ryT + 6'crr 
On = 

Similarly, 



{l-arf 



An+1,1 = /?n+i/?„+2 - crTa„+ia„+2 = (1 " ot)^", n > 
Substituting this and (5„ , 7„ into l|4.4(l , we get 

_ (l-^r)^ + (r? + ga)(g + r?T) 

^" - T. 'Ti^ ' ^ -J- 

(1 - err) 

We also compute 

(1 - o-r) 

Therefore recurrence (|2.19|l with initial values (|2.21|l gives 
xp2ix) = 1 - ar 1 + at)p3ix) + ^ '—^ p2{x) 

(1 — (TT)^ 

.(l-^^)' + (^ + ^^)(^tZ^(t + r)p,(.), 



(l-aT)3 



(?7 + 20a + ri(TT)t + 9 + 2riT + dar 



XPn{x) = (1 - (Tt)" (1 + (Jt)pn+l{x) H -3 p„(x) 

(1 - o-r) 

for n>2. 

After renormalizing the n-th polynomial by (1 — 0-r)("~^)("~i)/^, from Theorem 
14. II we get the following. 
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Proposition 4.3 (free quadratic harnesses). Suppose (Xt) is a quadratic harness 
with parameters such that (1 — crr)^ + (?? + Ocf){0 + rjT) > 0, < fr < 1, and 
q ~ —ar. If for t > random variable Xt has all moments and infinite support, 
then it has orthogonal martingale polynomials given by the three step recurrences 
n~Tl\ ), and 

xp2{x;t) = {l + at)p3{x;t) + ^ '—^ -2 P2ix;t) 

(1-ctt) 

+ ^ ^3 ^ t + r Pi , 

(1 — my 

, . . , , , {7] + 2ea + iqaT)t + e + 2r]T + eaT , , 

XPn{x;t) = [1 + C7t)pn+l[x;t) ^ "2 Pn\X]t) 

(1 - ctt) 

{l-aTf + {rj + ea){e + r,T) ^^ , / ,^ ^ o 

+ 77 {t + T)pn-i{x;t), n>2. 

(i — arp 

Remark 4.3. The recurrence in Proposition 14. 31 is a finite perturbation of the con- 
stant coefficient recurrence which was analyzed by many authors, see (29) and the 
references therein, see also |27| . 

Remark 4.4. In [2] we show that the free bi-Poisson process is associated with the 
generalized free convolution studied in 0. It is interesting to ask if an analogous 
situation occurs for the general free harnesses of Proposition 14.31 for recent exten- 
sions of generalized free convolutions to perturbations of higher-order terms, see 

m 

4.2. Classical quadratic harnesses. The classical quadratic harnesses have pa- 
rameter q = l — 2^/aT. The adjective "classical" is motivated by the fact that when 
err = 0, quadratic harnesses with g = 1 are related to classical probability, see [T^ 
Section 4.2]. 

Proposition 4.4 (Classical quadratic harnesses). Suppose {Xt) is a quadratic har- 
ness with parameters such that (1 — crr)^ + {rj -\- 0a){6 -f ryr) > 0, < err < 1, and 
q ~ \ — 2^JaT. If for t > random variable Xt has all moments and infinite support, 
then it has orthogonal martingale polynomials given by the three step recurrences 
17:171 ), 1777^ , and 

(4.16) xpn{x;t) = {l + {n-l)p + n(Tt)pn+i{x;t) + {-frit + 5n)Pn{x]t) 

(4.17) +ujl{{l + {n-2)p)t + T{n-^l))pn-i{x;t), n>2. 

where p = \faT , 

^- - (1 - rt^d X-^m. (2.. - m ^ (" " - "^^ 

(2n - 3)ryp -I- 2{n - ifr^p^ + (2(n - 1)^ - l)Ocfp^ , n > 1 
-1- {2n - i)ep + 2{n - ifOp^ + (2(n - 1)^ - l)ryrp) , n > 1 
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(4.20) n(l + {n-3)p) 



(1 - p)2(l + (2n - 2)p)(l + (2n - 4)p) 

n{n - 1)(1 + (n - 2)p)(l + (n - 3)p) 

^ (1 - p)4(l + (2n - 2)p)(l + (2n - 3)p)2(l + (2n - 4)p) 
X ( (1 + (n - 2)p) 6* + (n - l)r/T) ( (1 + (n - 2)p) ?] + (n - 1)6'ct) , n > 2. 

Proof. The assumptions of Theorem l4.1l are satisfied, so the orthogonal martingale 
polynomials exist, and we only need to solve the recurrences H4.2|l . (|4.3(l . H4.4|l . and 
renormalize the polynomials to simplify the final three step recurrence H4.17|l . 
From H4.2() we get 

a„ = (l-py-^n-l), n> 1, 
/3„ = (l-pr-2(i + („_2)p), n>2. 
Thus Ai^fc = f3k+i, k >0, and a calculation gives 
(4.21) A„,fc = (1 - p)2"+'=-3(l + (2n + A: - 3)p), n>2 

an+2f3n - Pn+2an = 2(1 - p)^"~"^, n > 1. 

Since g + ctt = (1 — p)^, equations 14.311 simplify to 

l + (2n-l)p 2(7 
4.22) 7„+i = — — — — p-7« + ■7-—Z — — p-On 
l + (2n+l)p l + (2n+l)p 

(1 + (2n - l)p + 2n'^p'^)T] + (2n + 1 + (2n'^ - l)p)9a 

+ (l-p)2(l + (2n+l)p) ' 

A l + (2n-l)p ^ 2r 

4.23 Sn+i = — — — — p-*^" + TTT^ — rrr^" 

l + (2n+l)p l + (2n+l)p 

(1 + {2n - l)p + 2n2p2)6> + (2n + 1 + (271^ - l)p)r]T 
+ (l-p)2(l + (2n+l)p) ■ 

Formula H4.10|l shows that H4.18|l and H4.19|l hold true for n — 1. Assuming that 
(|4.18|l and H4.19|l hold true for some n > 1, from equations H4.22|l and (|4.23l) a 
computer-assisted calculation shows that the formulas hold true for n + 1 as well. 
We now show that the solution of (|4.4(l is 

(4.24) CO,,- "(l + (n-3)p) 



(1 - p)2"-2(l + (2n - 2)p)(l + (2n - 4)p) 

n{n - 1)(1 + (n - 2)p)(l + (n - 3)p) 

(1 - p)2"(l + (2n - 2)p){l + {2n - 3)p)2(l + (2n - 4)p) 
X {{l + {n-2)p)9+{n-l)j]T){{l + {n-2)p)T]+{n-l)ea), n>l. 

Indeed, a calculation shows that the formula holds true for n — 1,2. Suppose (|4.24|) 
holds for some n > 2. For k = 1, formula (|4.21|) holds also for n = 1; thus for n > 2 
recurrence (|4.4f) simplifies to 

1 + (2n - 4)p 1 + T7,2 + aSl + Oj,, + 7?<5„ 2p(5„7„ 

(l-p)2(l + 2np)'^" (l-p)2"(l + 2np) 

Using I4.22|) . H4.23|l and 14.24II . a computer assisted calculation verifies that (|4.24|) 
holds true for n + 1 . 
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Renormalizing the n-th polynomial in (|2.19|l by the factor (1 — p)("~2)(n-i)/2^ 
we get (I4.17|l with w„ replaced by a;°. □ 

4.3. Orthogonal martingale polynomials when err = 0. We use the standard 
g-notation 

„n,-l 



= l + q-\ h 9" 



],! - [l]j2],...[n],, 



with the usual conventions [0]q = 0, [0]^! = 1. In this notation, Remark |4.1l gives 
an ^ [n- l]q, and /3„ = 1. 

Passing to the time-inverse (tXi^t) if necessary, without loss of generality we 
may assume that (7 = 0. In this case, recurrences in Theorem 14.11 have explicit 
solutions and orthogonal martingale polynomials are monic. 

Theorem 4.5. Suppose (Xt) is a quadratic harness with covariance ^2. b]) and 
parameters such that a = 0, —1 < q < 1, and 1 + [n]qr]9 + [n.]^T?7^ > for all n. If 
for each t > random variable Xt has infinite support, then the monic orthogonal 
martingale polynomials pn{x;t) are given by the recurrence 

(4.25) xpn{x;t) = pn+i{x;t) + {tjt + 6 + {[n]q + [n - l]q)r]T) [n]qPn{x;t) 

+ {t + T[n - l]q) (1 + [n - l]qri0 + [n - 1]^ V) [n],p„_i(a;; t), n>l, 
with the initial condition po = l,pi(x) = x. 

Remark 4.5. It is easy to give simple sufficient conditions on the parameters which 
imply that ojn := 1 + [n — l]qrj9 + [n — Ij^rr]^ > for all n > 1. Suppose that 
-1 < g < 1. 

(i) If 776* > then, trivially, ojn > 0. 

(ii) If r = 0, then a;„ > provided 1 + rjO > ma.x{q, 0}. Indeed, for n > 2 we 
have 

(4.26) l + g<N, <1 if-l<(/<0, 

(4.27) 1 <[n]q <l/{l-q) if < g < 1. 

For 77^ < we get w„ > from the right hand sides of these equations. 

(iii) if — 1 < (7 < 1 and 9^ < 4r then a;„ > 0. Indeed, then the quadratic 
function f{x) = 1 + xrjd + x'^rrj^ is non-negative. 

(iv) If r > 0, 776' < 0, \+T]9+Trj^ > 0, and 9^ > 4r, then the sufficient condition 
is 



(4.28) l + ,„„(,,0)>M±/fI3!. 

Indeed, inequality l + ri9 + Trj'^ > implies UJ2 > 0. Considering separately 
the cases 9 < and 9 > 0, the larger root of 1 + xri9 + x^rrf' = is 

'^'"^2rH~'^'" - If holds then the left hand sides of inequalities (g^ni, 

(|07|) imply uJn> for n > 3. 

It might be interesting to point out the explicit recurrence for the case r = 0. 

Corollary 4.6. Suppose {Xt) is a quadratic harness with covariance 1^2. b]) and 
parameters such that t = 0, —1 < q < 1, and 1 -|- [n]qr]9 + [n]qa9'^ > for all 
n. Suppose that for t > random variable Xt has infinite support. Then as the 
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orthogonal martingale polynomials for (Xt) we can take polynomials pn{x;t) given 
by the recurrence 

XPn{x]t) = {1 + ta[n]q)pn+lix]t) 

+ {0 + tr]+{[n]q + [n- l]q)0at)pn{x; t) +t{l + [n - l]qije + [n - l]^cr6|2) [n],p„_i (x; t). 

Proof of Theorem \4-5\ We use the notation from H2.2U|I . The integrabihty assump- 
tion of Theorem 14.11 is fulfilled by Theorem 12.51 and the system of recurrences in 
Theorem 14 . II simplifies . Since a„ = [n— 1]^, /3„ = 1, the first equation of recurrence 
(|4.3|1 becomes 

7«+i = '7 + InQ 

which gives 

(4.29) 7„ = r][n]g. 

The other recurrences are solved inductively. Suppose that 

(4.30) Sn-i = (6 + rjTi[n - 1], + [n - 2]g))[n - l]g , 

(4.31) LUn = {1 + [n ~ l]qf]e + [n - l]lTrj^)[n]q , 

hold. The initial conditions say that both formulas are satisfied for n — I, and the 
second one holds true also for n = 2. 

We use (|4.29|l and H4.3()|l to compute (5„ from the second equation in 14.3|l as 
follows 

Sn = + qSn-i + qi[n]q - - 2]g)T7„_i + r]T[n - l]q + riT[n]g = 0{1 + q[n - l]q) 

+7?r(g([n - 1], + [n- 2],)[n - 1], + q{[n]g - [n - 2],)[n - 1], + [n - 1], + [n],) 

= 6[n]q + i]T{q{[n - 1]^ + [n]q)[n - 1], + [n - 1], + [n]g) 

= 9[n]g + riT{[n - 1], + [n],)(l + q[n - l]q) = 9[n]q + T]T{[n - l]q + [n]q)[n]q 
Finally, (|4.4|l determines cun+i as 

UJn+l = qUJn + InT + QlrAi + InO + SnV + 1 - JnSn- 

Using the inductive assumption ()4.31|l , and already established formulas H4.29|l and 
ijOljl . we get 

uj„,+i = (1 + q[n]q) + (1 + q[n - l]q)[n]qr]0 + q[n - l]l[n]qTr]'^ + T7f[n]l 

+?7(1 + q[n]q)5n - ?/Ng(5„ ^ [u + l]q + [nYqTlO + q[n - lYq[n]qTTf- + TTf[nfq + Tjq'^Sn 

^[n+ l]q + [n+ l]q[n]qr^0 + Tr^^{q[n - l]q[n]q([n - 1], + q^-') + [n]^ + ^"[n]^). 
The coefficient at rrj^ simplifies to 

q[n - l]q[nfq + [nfq + = [n + l]q[nfq. 

Thus uj„+i = (1 + [n]qV0 + Ng + I]?- 

Formula H4.25(l comes from H2.19|l after substituting 14.2914.31(1 into 12.20(1 . 

□ 
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4.4. Operator solutions. In this section we re-derive the recurrences for some 
special orthogonal martingale polynomials from Section lT!^ bv an operator approach 
which is related to Lie algebra techniques. This method has a more ad hoc character 
so we concentrate on two relatively simple cases only. 

In the operator approach, we go back directly to Theorem 12. 31 We re-interpret 
the matrices X, y as the linear operators acting on the formal power series in an 
auxiliary variable z. This identifies the martingale polynomial pn{x;t) with z", 
where z is an auxiliary variable; similar technique appeared in umbral calculus j28l 
Ch. 1], in orthogonal polynomials ^Oli and in Segal-Bargmann representation |26j . 

We seek the solutions of in terms of the (/-differentiation operator 

9 [z) q=l, 

and the multiplication operator 

Z(.9)(z) = zg{z), 

treated as the linear operators on formal series g{z) in variable z. Table^lists the 
g-commutators of the combinations of these two operators that we need here and 
in Section H31 

The requirement on the 0-th column of reduces to the requirement that on 
the unit constant function Xl = 0, yl = z. When cr = 0, we have a„(t) = 1 so we 
are looking for the operators that satisfy 

(4.32) {tx + y)z" = + lower order terms. 





Z 


ZD, 


ZD2 


Z^D, 


D, 


1 


D, 




(l + g)ZD, -g(l-g)Z2D? 


ZD, 


z 


(l-g)ZD,+g(l-g)Z2D2 




Z^D, 


ZDl 


(l + g)ZD,-g(l-g)Z2D2 


ZD? 






Z^D, 


Z^ 









Table 1. q-commutators [A,B]q := AB - qBA. 



Example 4.7 {q-Meixnei processes). The g-Meixner process is a quadratic harness 
with parameters a = rj = 0, see jl2) . Inspecting Tabled we verify that 



and 



x = D, 



y = Z(l + ^D, + rD^) 



solve when a = r/ = 0. Indeed, 

[X, y], = [D„ Z], + 0[D„ ZD,], + r[D„ ZD% = I + ^D, + rD^. 
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A calculation shows now that 

{tx + y)z" = + e^gz" + {t + T[n - l]g)[n]gz"-\ 

which by the identification of z" with Pn{x; t) implies that the corresponding mar- 
tingale polynomials satisfy the three step recurrence 

(4.33) xpn{x;t) = pn+i{x;t) + d[n]qPn{x;t) + (t + r[n - l],)[n],p„_i(a;; i), n > 1. 

Of course, this is a special case of (|4.25|l corresponding to 77 = 0. Feinsilver |16[ Sec- 
tion 3.4] gives a reparametrization of recurrence (|4.33|l . and considers g-commutator 
relations [x, y]q = h as well as their realizations via operators D,, Z without consid- 
ering Anshelevich [3 Remark 6] discusses a reparametrization of recurrence 
(|4.33|l in relation to free Sheffer systems. In ^21 recurrence (|4.33|) is used as the 
first step in the construction of the g-Meixner Markov processes. 

We now use the same method to derive the recurrence for the martingale poly- 
nomials of the bi-Poisson process. 

Example 4.8 (bi-Poisson process). The bi-Poisson process is a quadratic harness 
with parameters a = t = 0. Inspecting Tabled we verify that 

x = D, + ryZ(D, + 0D^) 

and 

y-Z(l + 0D,) 
solve when a = t = 0. Indeed, 

[X, y], - [D„ Z], + ^[D,, ZD,], + r;[ZD„ Z], + r70[ZD„ ZD,], 

+r,0[ZDlZ]q + V0'[ZD' ZD,], = I + 0D, + ryZ 



+V9{1 - q)ZDq + rjOqil - q)Z^Dl + 7^6(1 + g)ZD, - ^Oqil - q)Z^Dl + -qOWq 
= I + 6ID, + 7?Z + 7?6'ZD, + r;6iZD, + r^e^ZDl 
= I + (0D, + riOZDq + ve'^ZD^q) + (?/Z + r;0ZD,) 

Operator tx -I- y satisfies the constraint {tx + y)l = z . A calculation shows that 

(tx + y)z" = + {9 + tr;)[n],z" + t(l + r/0[n - l],)[n],z""i, 

so the constraint (|4.32|) holds. By the identification of z" with pn{x;t), the corre- 
sponding martingale polynomials satisfy the three step recurrence 

(4.34) xpn{x; t) = Pn+i{x: t) + {9 + tr])[n]qPn{x; t) + t(l + rje[n - 1],) [n],p„_i (a;; t), 

n > 0, with p_i —Q^pq — 1. Of course, this is a special case of (I4.25|) corresponding 
to T = 0. But this recurrence was hard to guess without so in iH] it appears 

for (7 = only. 
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4.5. Dual g-commutation equation. Coherent states. Coherent states and 
the Segal-Bargmann representation are analytical techniques developed in mathe- 
matical physics j26) . The full Segal-Bargmann isomorphism is known to fail even 
in the relatively simple case of g-Brownian motion with q < 0, compare (31) . But 
algebraic duality is available and useful in a more general setting. 
Let 

(4.35) Qt* ,,„(x, y) = At,,,„x2 -I- Bt.s^uyy- + Ct.s^uf + A,.,nX -i- Et.s 

be the quadratic form in the non-commuting variables X,y; for background, see e.g. 
|21[ page 7] . Note that this is a dual of the quadratic form H3.32II that appears in 
the proof of Theorem 12.31 The following relates H4.35|l to the dual version of the 
g-commutation equation Hl.l(l . 

Proposition 4.9. Let Xt be a linear function Xt ~ X + ty of non- commutative 
variables X, y, and suppose that the coefficients of the quadratic form Ql s u '^''^ 
given by i2.15\} and \S.2,^8. £9)) . Then the following statements are equivalent. 

(i) Operator identity 

(4.36) X?-Q*,^JX.,X„) 

holds for all s < t < u; 

(ii) The non- commutative variables X, y satisfy the equation 

(4.37) [x,y]q ^ ax^ + Tf + rix+ey+l. 
Proof. Since X^Xs = x^ + suy^ + s(xy + yx) + (m - s)yx, we have 

= (1 - At.s,u - Bt.s,u - Ct,,,„)x2 + {t^ - S^At,s.u - SuBt.s,u - U^Ct.s,u)V^ 

+ (< - sAt,s,u - sBt^s,u ~ uCt^s,«)(xy + yx) - (m - s)St^^,„yx 

- (Dt.s^u + Et,s,u)y- - [sDt^s.u + uEt,s,u)V ~ Ft,s,J- 

Relations (|2.13|l . (|3.23|) . H3.24|l are equivalent to 

X? ~ Qt%,„(X„X„) = -<7Ft,s,u>^^ - TFt,s,uf 

= Ft,,,„([x,y], - ax2 - rf - ryX - 0y - I) = 0. 

□ 

4.5.1. Coherent states of q-Meixner process. Continuing Example 14. 71 the analogue 
of the coherent state in physics is the generating function 

(4.38) ipt[z, a;) = V Y-—pn{x\ t), 

where p„ (x\ t) are the orthogonal martingale polynomials given by recurrence (|4.33|l . 
Proposition 4.10. The operator 

(4.39) Xt = {l + eZ+TZ^)Dg + tZ 

satisfies i4.cl&\) with the quadratic form Ql g.u g'iven by l\2.15}l . 1^3.2513. 2y\) where 
a — rj — Q. Moreover, 

(4.40) x'pt[z,x) = {Xt(ft){z,x). 
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Proof. Inspecting Tabled it is easy to verify that X = {1 + 9Z + rZ^)Dg and y — Z 
satisiy H4.37|l with cr = r; = 0. Indeed, 

[x,y], - [D,,Z], + e[ZD,,Z], + r[Z'D,,Z], 

= \ + eZ + TZ^ = l + 6iy + Ty^ 

Therefore, Proposition l4.9l imphes H4.36|l . The algebraic identity 14.4UII follows from 
recurrence (|4.33|l by the following calculation. 

°° z" 

xipt{z,x) -— -xp„ {x; t) 
= rn iPn+i{x;t) +9[n]pn{x;t) +T[n][n- l]p„_i(x;t) + t[n]pn-i{x;t)) = 

n— 1 '- n— 1 '- 

= (D, + 6iZDg + rZ^Dg + tZ)ipt{z, x). 

□ 

Formula H4.40|l is implicit in the usual derivation of the product formula for 
(pt{z,x), compare [J. When the parameters t,0 vanish, it appears in |31l Section 

31 in the context of analyzing ground states for the g-deformed Gauss distribution. 
Ref. 12] gives a more general analytical scheme which coincides with H4.39|l and 
(I4.40|l when t = 0,q — 1; however, it advocates the normalization by the L2-norm 
of the polynomials, which does not fit all the cases we are interested in. 

In we defined the g-Meixner process as a Markov process with the initial state 
Xq = and with the transition probabilities Ps^t{x,dy) determined as the unique 
probability measure orthogonalizing the polynomials Qn in variable y which are 
given by the three step recurrence 

(4.41) vQn{v\x) 

= Q7i+i{y\x,t,s) + {0[n]q + xq'^)Qniy\x,t,s) 
+ {t - sq""-^ + T[n - l]g)[n]gQn^i{y\x, t, s). 

In that paper, we showed that this Markov process is well defined, that it has 
orthogonal martingale polynomials Pn{x\ t) given by recurrence H4.33|l . and we used 
this to prove that {Xt) is a quadratic harness with parameters a — rj — Q. 

Proposition 14. l()l simplifies the verification of the quadratic harness condition in 
Proposition 3.4], condensing more than three pages of proof into one page. 

Proposition 4.11 (jEj). 1} \q\ < 1 and the polynomials Pn{x;t) defined by i4.3S[ ) 
are orthogonal martingale polynomials for a Markov process (Xt), then (Xt) is a 
quadratic harness with parameters cr = r] = 0. 
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Proof. For simplicity, we consider only |g| < 1, as in this case H4.Hr>ll implies that 
\Xt\ < Ct has bounded support. Consider the generating function H4.r>8|) . Since 
(Xt) is Markov, condition H2.3|l is equivalent to 

(4.42) E{ip,{zi,X,)Xt^u{z,X,,)) 

= at.s,u^ {^sizi, Xs)Xs(psiz, Xs)) + bt.s,u^ {ipsizi, Xs)Xu(puiz, Xu)) 

holding for all zi, z in a, neighborhood of (or just as an identity in formal power 
series in variables zi, z). 

We now use the martingale polynomial property which implies that 

(4.43) E[Mz,Xt)\X,]^^siz,X,), 

and we use (|4.40|l to represent the process through the operator H4.39(l as 

Xtiptiz,Xt)^XtiMz,Xt)). 

This gives 

E {^s{zi,Xs)Xt^u{z, Xu)) = E {^s{zi,Xs)Xt^t{z, Xt j) 
= '>^t^{'Ps{zi,Xs)(ps{z,Xs)) = XtGs{zi,z), 

where 

(ziz)' 

n=0 



iz zl" 

Gs{zi,z) = E{^s{zi,Xs)^s{z,X,)) = t47^E(p„(X,;s)2) 



n=0 ^ J* fc=l 



and Xt acts on Gs{zi,z) as a series in variable z. Therefore, equation ()4.42|l is 
equivalent to 

XtGs(zi,z) = at^s,uXsGs{zi,z) + 6t^s,„X„Gs(zi, z), 

and follows from the (trivial) operator identity Xt — at,s,uXs + bt,s,uXu- 
Similarly, condition H2.11II is equivalent to 

{A.U)E{<Ps{zi,Xs)X^<Pu{z2,Xu))^At,s.uE{ips{zi,Xs)Xl^s{z2,Xs)) 

+Bt,s,u^ {iPs{zi,Xs)XsXu(Pu{z2,Xu)) + C(,s,„E [ips{zi,Xs)Xlipu{z2,Xu)) 
+A,s,mE {ips{zi,Xs)Xs(fis{z2,Xs)) + £^t,s,«E {ips{zi, Xs)Xu(fiu{z2, Xu)) 
+Ft.s,nM' {'Ps{zi:Xs)(Ps{z2,Xs)) ■ 

Notice that for s < u we have 

E [ips[zi,Xs)XsXu'Pu[z, Xu)) = E {(fsizi, Xs)XsXuipu{z, Xu)) 
= XuEiifisizi, Xs)Xsipsiz, Xs)) XuXsEiips{zi, Xs)iPsiz, Xs)) = XuXsGs{zi,z). 

Therefore, equation H4.44|l follows from the operator identity H4.36|l which in ex- 
panded form says 

X( — At^s,uX^ + Bt^s,uXuXs + Ct^s.uX^ + Dt^s,uXs + Et,s,uXu + Ft^s,u\, 

applied to Gs{zi, z) treated as a series in variable z. □ 
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4.5.2. Coherent states of the bi-Poisson process. We now repeat the methods of 
Seetion l4.5.1l to derive new resuhs about the bi-Poisson process from Example 14.81 

Proposition 4.12. The operator 

(4.45) Xt = {l + {e + i]t)Z+ tr]0Z^)Dq + tZ, 

satisfies i4.36}) for the quadratic form Qt s u with the coefficients \2.15]) . J^H. 2^3. W^) 
such that a ~ T ~ Q. Moreover, if (pt{z, x) is the generating function i4-'^^ of the 
orthogonal martingale polynomials Pn{x;t) given by recurrence 

then 

(4.46) xiptiz,x) = {Xt(pt)iz,x). 
Proof. Inspecting Tabled we verify that 

X= (1 + 6IZ)D, 

and 

y = Z + 77(Z + 0Z2)D, 
solve (|4.37|l with cr = r = 0. Indeed, 

[X, y], = [D„ Z], + d[ZD„Z], + 77[D„ ZD,], + 7/0[ZD,, ZD,], + ry0[D„ Z^D,], 
+V0^[ZDg, Z^D,], = \ + 0Z + 7^Dg + f]6{{l - g)ZD, + g(l - q)Z''Dl) 
+Tje{{l + q)ZDg - q{l - q)Z^Dl) + rjO^Z^Dg = I + 0Z + 77D, + ry^ZD, 

+ry0ZD, + rje^Z'^Dq ^ \ + r]{\ + 0Z)D, + 0(Z + r;(Z0Z2)D,) ^ I + 77X + 0y. 

Therefore, Proposition 14. 91 implies (|4.36|1 . 

We now derive (|4.46|l from (|4.34|l by the following calculation. 

°° z" 

XLpt{z,x) = V --—-xpn(x;t) 
^ pp {Pn+i{x;t) + {6 + tr])[n]qPn(x;t) + t{l + r]e[n - l],)N,p„-i(a;; 0) = 

n=0 



^^lYl r I 11 , Pn+i{x; t) + {e + iry)zD, ^ -— p„(x; t) 

„=0 ^'^"^ -^^9- n=0 

r _ 11 | P»-i(a:;^) +i?76'z^D, ^ r _ 11 | P»~i(^;^) 
= ((1 + {9 + Tjt)Z + t7j9Z^)Dg + tz) ipt{z, x). 



□ 



Next we give the bi-Poisson version of Proposition 14.111 



Proposition 4.13. Fix — 1 < q < 1. If polynomials Pn{x;t) given by i4.34}) are 

orthogonal martingale polynomials for a Markov process (Xt), then (Xt) is a bi- 
Poisson process, i.e., ^2.3}) and ^2.14^ hold with a = t = 0, and l + rjO > max{g, 0}. 
Moreover, such Markov process (Xt) is determined uniquely. 



QUADRATIC HARNESSES 



33 



Proof. Property (|2.6|1 follows from the explicit form of the polynomials po,Pi,P2- 
Orthogonality gives E{Xt) = E{pi{Xt; t)po{Xt;t)) = and E{X?) = E((p2(^t; t) + 
{d+rit)Xt+t)pQ{Xt)) =t. Martingale polynomial property then implies E(XsX4) — 
E(X,E[pi(Xt;t)|^<,]) =E(X2) 

If —1 < q < 1, from (|4.34() it follows that|Xt| < Ct has bomided support; if 
q = 1 then it is not hard to identify the distribution of \Xt\ for example, from 
Chapter 4] or pages 175-181], verifying that \Xt\ has a finite exponential 
moment. Thus the moment problem has unique solution and the process (Xt) 
is determined uniquely. Moreover, polynomials are dense in L2{Xs, X^^), see |15l 
Theorem 3.1.18]. 

Consider the generating function H4.38|l with polynomials Pn{x;t) satisfying 
(j4.34|l . Notice that property 14.43|l follows again from the martingale polynomial 
condition. By ProDOsition l4.12l with Xt defined by l|4.45(l we have the representation 

Xtipt{z,Xt)^XtiMz,Xt)). 

Since (Xt) is Markov and polynomials are dense in L2{Xs, Xu), condition (|2.3() 
is again equivalent to l|4.42|l . which we can interpret as the identity between the 
formal power series in variables zi,z. The latter follows again from the (trivial) 
operator identity Xj = at_s_„Xs + bt^s,uXu, applied to 

G,(zi,z) =E(^,(zi,X,)(^,(z,X,)) = J2 TinT^(P"(^-;^)')' 

n=0 

treated as the formal power series in variable z. 

Similarly, condition 12.1111 is equivalent to (|4.44|) . and again 

E {ips{zi, Xs)XsXu(pu{z, Xu)) ^ XuXsGs{zi,z). 

Therefore, equation (|4.44|l follows from the operator identity (|4.36(l applied to 
Gsizi, z) treated as the formal power series in variable z. 

Since the third coefficient in H4.34|l is nonnegative for all n, we get 1 + rjO > 
maxiq. 01. compare Remark l4.5l □ 
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